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Abstract 

In this work we study an SO(10) GUT model with minimum Higgs repre- 
sentations belonging only to the 210 and 16 dimensional representations of 
SO(10). We add a singlet fermion S in addition to the usual 16 dimensional 
representation containing quarks and leptons. There are no Higgs bi-doublets 
and so charged fermion masses come from one-loop corrections. Consequently 
all the fermion masses, Dirac and Majorana, are of the see-saw type. We min- 
imize the Higgs potential and show how the left-right symmetry is broken in 
our model where it is assumed that a D-parity odd Higgs field gets a vacuum 
expectation value at the grand unification scale. From the renormalization 
' group equations we infer that in our model unification happens at 10 15 GeV 

■ and left-right symmetry can be extended up to some values just above 10 11 

GeV. The Yukawa sector of our model is completely different from most of the 
standard grand unified theories and we explicitly show how the Yukawa sector 
will look like in the different phases and briefly comment on the running of the 
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top quark mass. We end with a brief analysis of lepton number asymmetry 
generated from the interactions in our model. 
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1 Introduction 



The SO (10) grand unified theory has several interesting features [1-5]. It can accom- 
modate left-right symmetry as one of the intermediate symmetry and hence provides 
an explanation of parity violation [6]. (B — L) is a generator of the group SO(10) and 
hence lepton number violation takes place spontaneously. This explains the origin 
of lepton number violation and neutrino Majorana mass naturally. The smallness 
of neutrino masses is assured by the see-saw mechanism, so that by keeping the 
scale of (B — L) violation high the smallness of the neutrino mass is guaranteed. 
Quarks and leptons are treated equally in SO(10) GUT. Gauge coupling unification 
is consistent with all low energy results. 

On the one hand there are many attractive features of the SO(10) GUT, on the 
other the predictability becomes low. Depending on the symmetry breaking pattern 
and Higgs scalar contents, the model can have widely differing predictions. Attempts 
have been made to construct a minimal model. In one approach the minimal model 
is constructed with minimum numbers of parameters, while in the other approach 
minimum numbers of Higgs scalars are included in the model. There are also models 
without any intention of minimality or simplicity, where the main aim is to explain 
all experiments and have maximum predictability. 

In the present article we shall study an SO(10) GUT, which has the minimum 
dimensions of the Higgs scalars. In any SO (10) GUT the minimal number of Higgs 
scalar includes an SU(2) L symmetry breaking Higgs scalar, which can give masses to 
the fermions and some scalars that can break the (B — L) symmetry along with the 
SU(2) R symmetry and can give Majorana masses to the neutrinos. In addition, there 
is one Higgs scalar which breaks the group SO(10) GUT. One conventional model 
includes a Higgs bi-doublet (a 10-plet of SO(10) Higgs scalar, which is doublet under 
both the SU(2) L and SU(2) R groups) and both SU(2) L and SU(2) J? Higgs triplets. 
In such models with triplet Higgs scalars neutrinos acquire masses at the tree level. 
The SO(10) representation that contains this Higgs scalar is of 126 dimensions. In 
another version of the model, one breaks the left-right symmetry and the (B — L) 
symmetry by doublets of SU(2) L and SU(2) R groups. This Higgs field belongs to a 
16-plet representation of SO(10). For symmetry breaking and giving fermion masses 
another 10-plet Higgs scalar is introduced, which contains the bi-doublet Higgs and 
hence gives tree level masses to the fermions and the neutrinos acquire only Dirac 
masses of the order of other fermion masses. There is no Majorana mass term for 
the neutrinos and hence see-saw mechanism is not possible. However, there exist 
effective higher-dimensional operators, which can give correct Majorana masses to 
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both the left-handed and right-handed neutrinos. 

Recently it has been pointed out that it is possible to consider an SO(10) GUT, 
which does not have any Higgs bi-doublet scalar belonging to a 10-plet of SO(10) [7]. 
The Higgs scalar that breaks left-right symmetry and the (B — L) symmetry belongs 
to a 16-plet of Higgs scalar. Since tree level fermion masses are not allowed without 
the bi-doublet scalar, all fermion masses come from higher-dimensional operators 
in the see-saw form. In supersymmetric theories the non-renormalizablity theorem 
does not allow radiative generation of such higher dimensional operators. We shall 
then restrict ourselves to only non-supersymmetric SO(10) GUT. The source of see- 
saw suppression for the fermion number violating Majorana mass terms are different 
from the source of see-saw suppression for the fermion number conserving Dirac mass 
terms, which maintain the large hierarchy between the charged fermion masses and 
the neutrino masses. In this article we shall study some aspects of this model in 
detail. 

In the next section we shall discuss the model. In sec. |3]we shall present details 
of the scalar potential minimization and the allowed symmetry breaking pattern and 
in sec. H] the generation of fermion masses is discussed. We shall then study the 
renormalization group equation for this model with the specific choice of the Higgs 
scalars. In sec. we shall study the gauge coupling unification and in sec. EH the 
Yukawa coupling evolution for the different fermions. Since all fermion masses have 
the same see-saw origin, the perturbative unification becomes an important question 
in these models. Some of the Yukawa couplings could become large, although the 
effective fermion masses still remains small. In sec. [7| leptogenesis in our model is 
discussed and in the last section we summarize our results. 



2 The Model 

The starting point for any SO(10) GUT is the choice of the symmetry breaking 
pattern. There exists many chains of symmetry breaking pattern, which are all 
consistent with our present knowledge. So the particular choice of a symmetry 
breaking pattern defines a specific model. We shall consider a symmetry breaking 
pattern, which requires a minimum number of Higgs scalars, given by: 

SO(10) ^ SU(4) C x SU(2) L x SU{2) R [G 422 ] 

SU(3) C x SU(2) L x S\J(2) R x U(1) B _ £ [G 3221 ] 
^ SU(3) C x SU(2) L x U(l) y [G321] 

^ SU(3) c xU(l) em , (1) 
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where the Higgs fields responsible for the symmetry breakings, <3> 4 , xr an d Xl 
are explicitly shown in the above equation. Both $£> and $ 4 are contained in $, 
which transforms as the 210 dimensional representation of SO(10). The 210-plet 
decomposes under the Pati-Salam subgroup G 42 2 of SO(10) as, 

$ = 210 = (1,1,1) + (6,2,2) + (15,3,1) + (15,1,3) 

+ (15, 1,1) + (10, 2, 2) + (10,2, 2). (2) 

In the above decomposition $d corresponds to (1,1,1) and $ 4 corresponds to 
(15,1,1)- 

The Higgs fields xl an d x*r belong to the 16 dimensional spinor representation 
T of SO(10) and the other fields xr an d Xl belong to the conjugate representation 
16, called 1^, of SO (10). The group transformation properties of the x fields under 
G3221 and G 42 2 are as follows: 

Xl = (1,2,1,-1)C(4,2,1)C16 

X* R = (1,1,2,1)C(4,1,2)C16_ 

Xr = (1,1,2,-1) C (4,1,2) C_16 

Xl = (1,2,1,1)C(4,2,1)C16 

In the above equations (x, y,z,w), and (x, y,z) denote the group transformation 
properties of the Higgs fields under G3221 and G 4 22- 

At this stage we shall digress to discuss one important feature of the left-right 
symmetric models, namely the question of parity V. The discrete Z2 symmetry, 
that interchanges the two SU(2) subgroups of the Lorentz group 0(3,1), is called 
the parity. This parity can be identified as the discrete Z 2 symmetry operator 
that interchanges the groups SU(2) L and 811(2)^ of the left-right symmetric model, 
which implies that under parity Wl < — ► Wr. This definition extends to scalars 
also. That is, an SU(2) L doublet scalar field xl will transform to an SU(2) R doublet 
scalar field xr under the operation of parity xl i — > Xr- m the conventional left- 
right symmetric models, the parity is spontaneously broken along with the group 
SU(2) R . In other words, when the left-right symmetric group SU(2) L x SU(2) R is 
spontaneously broken, parity is also spontaneously broken. 

There is another possibility of breaking parity spontaneously without breaking 
the left-right symmetric group. Since the scalar fields transform trivially under the 
Lorentz group, the VEV of a parity odd field can break parity spontaneously without 
breaking the left-right symmetry. Unlike the conventional case, now the parity acting 
on the fermions and vector bosons is not spontaneously broken. To distinguish these 
two cases, this second type of parity is called a -D-parity. Thus when D-parity is 
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broken, the left-handed and right-handed scalars can have different mass and VEV 
and hence the gauge coupling constants of SU(2) L and SU(2) H can also be different. 
In the present model D-parity plays a very crucial role, both for symmetry breaking 
as well as for fermion masses. It also plays some role in gauge coupling unification. 

The 210 representation of SO(10) is a totally antisymmetric tensor of rank 
four <&abcd and the singlet $d is the component $6789 m the notation, in which, 
a,b,c,d = 0,1, ..,5 are SO(6) indices and a,b,c,d = 6,7,8,9 are SO(4) indices. 
Thus under D-parity & D is odd ($_d — > — ®d) and consequently when it gets its 
vacuum expectation value (VEV) at the GUT scale, M U: it breaks the left-right 
parity of the theory. Due to this spontaneous breaking of the left-right parity at the 
M v scale we will have (xl) 7^ (xr) at a lower energy scale. This D-parity odd field 
is also required to give masses to the light neutrinos. 

Next we write down the fermions in our model and their group transformation 
properties. The left-handed quarks, leptons, anti-quarks and anti-leptons belong to 
a 16-plet representation of SO (10), which transform under G422 as: 



i = 1,2,3 is the generation index. The right-handed fermions and anti-fermions 
belong to the conjugate representation, 



In addition to the above mentioned conventional particles our model consists of an 
extra SO(10) gauge singlet fermion per generation: 



16 =(4, 2,1) + (4, 1,2), 



(4) 



i, iR = 16 = (4, 2,1) + (4, 1,2). 



(5) 



S iL = l = (1,1,1), 



(6) 



1 = 



1,2,3. 

Under G3221 the states (4, 2, 1) and (4, 1, 2) transform as 



(4.2.1) = (3, 2, 1,1) + (1,2, 1,-1), 

(4.1.2) = (3, 1,2,-1) + (1,1, 2,1), 



(7) 



(8) 



and as a result the fermions can be labelled as: 




(10) 
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and 

q R c = q c L = ; " I A 1.2. - (11) 



-R 




;i2) 



The generators of the left-right symmetry group G3221 are related to the electric 
charge of the particles by, 

Q = T 3L + T 3R + = T 3L + I , (13) 

where 

Y = T, R+ i ^±. (14) 

In the conventional left-right symmetric models there is one bi-doublet Higgs scalar 
= (1,2,2,0), which gives masses to quarks and charged leptons and a Dirac 
mass to the neutrinos through its couplings of the form ip L ip R (f). In addition, there 
are triplet Higgs scalars A^ = (1,3,1,-2) and A R = (1,1,3, — 2), which can 
give Majorana masses to the left-handed and right-handed neutrinos through the 
couplings ipL^ipL^^Lfi- In our present model all these Higgs scalars 0, Al and A^ 
are absent and hence there are no tree level fermion masses for the quarks and the 
leptons. After discussing the structure of the Higgs vacuum expectation values in 
this model, we shall come back to the question of fermion masses. 



3 Minimization of the scalar potential and left- 
right symmetry breaking 

In the conventional left-right symmetric models, the combinations of the Higgs fields, 
0, Al and A R ensures that for certain choices of parameters, A R can acquire a very 
large VEV compared to other fields breaking left-right symmetry at a large scale. 
It is clear that in the absence of the field 0, both the fields A^ would acquire 
equal VEVs. It has also been shown that in the absence of the field in a left-right 
symmetric model with only the doublet Higgs scalars xl and \R) the minimization 
of the potential would result in equal VEVs for both xl and XRi which would lead 
to inconsistency and parity will be conserved at low energy. This problem could be 
solved if parity is broken in these theories either explicitly or spontaneously. 

In the present model this problem does not occur. It was mentioned in the 
original version of the model that the D-parity odd singlet field <3>d = (1,1,1) 
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under G 42 2 contained in the 210 representation would allow (xr) (xl) and break 
left-right symmetry at some high scale. In this section we shall minimize the scalar 
potential and discuss the various possible solutions, which allows left-right symmetry 
breaking at some high scale. 

Let us consider the Higgs potential [7]: 

c s = m|$ 2 + a$ 3 + ^$ 4 + m2(rtr) + ^(rtr) 2 

+ ^[r 4 + (rt) 4 ] + M D ${r*r) + A^^r) . (15) 

The coupling ^fT^T) is the most important term that is required for the left-right 
breaking to take place at a higher scale compared to standard model symmetry 
breaking. D-parity is broken when $£> acquires a non- vanishing VEV, ($d) = rj, 
at the Mjj scale, since <J>d is odd under D-parity. $4 will get a non- vanishing VEV 
at the Mx scale. There will be many terms in Eq. (fTKJ) including ($4), but as we 
are analyzing the structure of Eq. (JT5j) in the G3221 phase and mainly interested on 
the VEVs of the x fields, we do not explicitly write down the terms including ($4). 
($4) has no important contribution in the expressions of the VEVs of the x fields. 

We now discuss the masses of the components of F and the VEVs. The scalar 
potential responsible for the masses of the fields xl and xr is given by, 

V = m 2 A\XR\ 2 + \Xl\ 2 ) + M D V (\x R \ 2 - \ X l\ 2 ) 

+ ^tV 2 (\Xr\ 2 + \Xl\ 2 ). (16) 

The masses of these fields are then given by, 

A = m r ~ m d V + -W V 2 , 

/4 = m 2 r + M D rj + A$ r . (17) 

If /^-parity is conserved, rj = and the masses of both xl and xr become equal. 
Since the VEV of $d breaks D-parity, it will be possible to fine tune parameters to 
obtain the mass of xl to be orders of magnitude smaller than the mass of xr- From 
phenomenological consideration we also require 

(Xl) = v L ~ n L ~ 100 GeV , 

(Xr) = v R ~ [x R ~ Mu^>v L . (18) 

We shall next check if this widely different VEVs for xl and xr is possible, vl 
breaks the electroweak symmetry, while vr breaks the left-right symmetry at a very 
high scale, close to the GUT scale. 
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We denote the VEVs of the fields xl and xr as: 



(Xl) 



v l 




, (Xr) 



, (Xr) 



Vr 




, (xl) 



(19) 



Instead of minimizing the potential, we shall first write down the potential in terms 
of the VEVs of the fields and then find the conditions satisfied by the VEVs. With 
the above VEVs we can write the Higgs potential in the G3221 phase as: 



V 



X 



4! 



+ [ml + X w (v 2 + ■ ■ -)](v R v R + v* L v L ) - ^( V R V R + v *l v l? 



X{ 



r (t7 4 + v 4 + w ,4 + + MDr){v * RVR _ v1vl) . 
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(20) 



The • • • symbols in the above equation stands for terms containing ($4). 

Setting vr = v* R and vl = v* L , which amounts to saying that there is no CP 
violation and all VEVs are considered to be real, the extremum conditions of V 
comes out to be: 



dV 

dv L 

dV 

dv R 



2vi 



2v 



R 



' " + v 2 R ) - X' r v 2 L - M DV 

ml + (X^rri 2 + ■■■)- y (^ + v%) - \' r v 2 R + M m 



0, 



0. 



The above equations imply, 

( dV\ ( dV 



dv 



vr 



dvi 



2v L v R [X' r (v 2 L - 4) + 2M D r]} = 



(21) 



(22) 



(23) 



Neglecting the trivial solution vl = v R = 0, the other interesting relation between 
vl and v R that comes out from the above equation is, 

2M D r] 



v r-v 1 



(24) 



Two things can be noted from the above equation. First as it was stated previously, 
in understanding the relation between vl and v R we do not require the VEV of $4. 
Secondly if rj has some value comparable to My and X' T is not too high, then it is 
apparent from Eq. (J2~4~j) that v R ^> Vl- If the energy scale where Xr an d X*r gets 
a non vanishing VEV be M R then we can say that M R 3> M\y where M\y ~ 100 
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GeV. Thus this model allows left-right symmetry breaking at a much higher scale 
compared to the standard model symmetry breaking scale. 

It is clear from the above discussions that this model works only if D-parity 
is broken spontaneously. In addition, severe fine tuning is required to obtain and 
maintain this solution. To make the masses of xl and xr. different /ii ^ fi R , a fine 
tuning is required. Then the next fine tuning is required to keep the VEV Vl to be 
orders of magnitude smaller than v R . This is the usual fine tuning required in all 
non-supersymmetric theories. We can write Eq. (f2~T| as 



dV 

dv L L 



2 i 2 i 2\ \/2 

+ v R ) - X r v L 



. (25) 



Since the VEV vl will be proportional to /i^, a fine tuning is performed to keep 
fiL ~ 100 GeV. The second fine tuning makes sure that the VEV v R does not 
destabilize the VEV of vl through radiative corrections. 

4 Fermion masses 

In the left-right symmetric theories the left-handed fermions are doublets under 
SU(2) L and the right-handed fermions are doublets under SU(2) R . Hence the fermion 
masses would require a bi-doublet Higgs scalar. Following our discussions at the 
end of sec. it is clear that in the present model there are no Yukawa couplings 
giving Dirac or Majorana masses to the quarks and leptons. In this model both 
the Majorana and the Dirac masses originate from dimension-5 effective operators, 
given by: 

Ox = -^-{WXL){qRX*R) , 2 = -^-{qExt)(QRXR) , 

1V1 D M D 

Os = -^-(hxDihxl) , = -L(£ rXr )(£ rXr ) , 

M M m m 

where Mp, M' D and Mm are some heavy mass scales in the theory. In general, the 
mass scales appearing in the operators which contribute to the Dirac masses (Md, 
M' D ) and the mass scales that appear in the operators contributing to the Majorana 
masses {Mm) will be different, since in one of them total fermion number is violated 
by 2 units. 
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When the Higgs scalars xl and xr acquire VEVs, the first two operators 0\ 
and O2 give the quark masses: 

,n, Wm , with m ,= ^. (27) 

Similarly the third and the fourth operators 3 and O4 contribute to the charged 
lepton and neutrino Dirac masses: 

mi l L £r, with m e = , (28) 

while the last two operators O5 and Oq contribute to the Majorana masses for the 
left-handed and right-handed neutrinos respectively: 

m ^ = lur~ and m ^ = Ivh ■ ( 29 ) 
Mm M m 

We shall now discuss some of the possible origin of these operators and their conse- 
quences. 

The see-saw masses of the neutrinos in theories with only doublet Higgs may 
arise from various cases as, some higher dimensional effective operators in a non 
super symmetric theory, from non-renormalizable gravitational interactions or from 
super symmetric extensions of models with doublet Higgs [10]. In the present case 
the see-saw masses of the neutrinos can be obtained in three different ways. They 
may be mediated by exchange of scalar fields or fermion fields or may be induced 
radiatively. As we shall argue now, the first two possibilities are not very attractive 
and hence we shall study the radiative mechanism in details. 

When the intermediate field is a scalar, it has to be a field which transforms 
as 16 x 16 and hence the field could be either a 10 or a 120 or a 126. If the scalar 
field transform as 120, the fermion mass matrix will be totally antisymmetric and 
hence phenomenologically unacceptable. If the scalar field £ transform as a 10 or 
a 126, its components will receive induced VEVs through its couplings XTT and 
Sr^T'''. Then we can eliminate the xl and xr i n the resulting theory and revert 
to the conventional theories with bi-doublet Higgs and triplet Higgs scalars A L R . 
So, we shall not discuss this possibility any further in the rest of the article. 

We shall now consider the possibility of intermediate heavy fermions generating 
the effective operators for the quark and lepton masses. For each of the operators we 
require two fermions, one left-handed and the other right-handed, both having same 
gauge transformation properties. For the Majorana mass terms generated by the 
last two operators a self-conjugate singlet fermion is sufficient. The singlet fermion 
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SiL, we already included in the present model, can give the Majorana masses to the 
left-handed and right-handed neutrinos. 

To generate the operator 1 , we need two fermions Ul and U R coupling to qlXL 
and qnX*R respectively. Both these fields should then transform similarly Ul,r = 
(3, 1, 1,4/3) C (15, 1, 1) C 45 or 210 and the Lagrangian must contain the couplings 

£i = aiU L q R x* R + biqZU R XL + muuU L U R + h.c. , (30) 

to give masses to the up-quarks by the operator 0\. The down quark masses are 
obtained by an effective operator 2 , which may be generated by adding the field 
Dl,r = (3,1,1,-2/3) C (6,1,1) C 10 or 126 or D LiR = (3,1,1,-2/3) C 
(10, 1, 1) C 120 and introducing the couplings in the Lagrangian: 

C 2 = a 2 D L q R x R + b 2 qZD R x L + m DD D L D R + h.c. . (31) 

The operators O3 and O4 may be obtained by introducing the fields Nl,r 
; 1. 1. 1.0; c (1,1,1) C 1 or 45 and El, r = (1,1,1,-2) C (10,1,1) C 120 
with the couplings 

C e = a 3 N2e R x R + b 3 hN RX L + m NN N2N R 

+ a 4 E L £ R x R + b A i L E R x* L + m EE E L E R + h.c. , (32) 

respectively. Then we may give masses to the up and the down quarks as well as 
to the charged leptons and the neutrinos if there are heave fermions transforming 
as 120 and 45. The singlet field Sn per generation is required to give Majorana 
masses to the neutrinos with its couplings, which we shall discuss later. 

We shall now come back to the present model, where the quark and lepton 
masses are generated radiatively. The fermion content of the model has been dis- 
cussed in sec. 121 The most general Yukawa couplings are then given by, 

C Y = fS L ^ + M S S L S L + h.c. . (33) 

In this expression generation indices have been suppressed. One loop diagram of 
Fig. ^ then generates effective operators 

4>l 4>r xl x* r c ipL V'r r r , (34) 

which are of the form 0\ and 3 and contributes to the down-quark and charged- 
lepton masses. On the other hand the one loop diagram of Fig. El generates effective 
terms: 

ipL iPr Xr Xl C ip L ipR r t r f , (35) 
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C 11 

u S S 11 



L 



R L R 

Figure 1: One loop diagram contributing to the fermion masses. 

which are of the form of the operators Oi and O4 and contributes to the masses of 
the up-quarks and the Dirac masses of the neutrinos. 

The up-quark, down-quark and charged-lepton masses can now be estimated 
from Fig. Q]and Fig. [21 to be: 

M u = (36) 
8vr 2 M x 

M df = (37) 

Here Mx = .1 ly / A Is or Ms, depending on whether Mr or Ms is larger and tul = fv^, 
and rriR = Jvr. We thus obtain different up and down quark masses and on the 
other hand b — r unification. The other mass relations in the down-quark sector 
and the charged-lepton mass relations could come from higher order terms, since 
the remaining matrix elements are of the order of 10~ 3 to compared to the 
33-element [8,9]. For example, operators of the form 

fe^rr$ 4 and J~ L i) R r f r f $ 4 , (38) 

contribute differently to the down-quark and charged-lepton masses, since the effec- 
tive VEV transform as rr$ 4 and r' $ 4 , which behaves as the field (15, 2, 2) C 126 
and hence can solve the fermion mass problem in GUTs, a la Georgi-Jarlskog mech- 
anism. 

The neutrino masses come from the couplings of the neutrinos with the singlet 
fermions Sn, given by Eq. In the basis ( v L v c L ) the tree level neutrino 
mass matrix becomes: 

/ m L \ 
M v = m R \, (39) 
\ m L m R M s I 
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Figure 2: One loop diagram contributing to the fermion masses. 



which gives two heavy states, which are mostly Sl and v c j J . In the limit M$ 3> mjj, 
the two heavy mass eigenvalues are Ms and m\jMs- On the other hand, when 
Ms <C ttir, the two heavy states are almost degenerate with eigenvalues ±7Br 
with a mass splitting of about M5. The latter case may be more interesting for 
leptogenesis, which we shall discuss at the end. 

The lightest state vl remains massless at the tree level. However, if we include 
the effect of D-parity violation, this problem could be solved. We thus continue our 
discussion taking D-parity violation into consideration. The effective operator: 



v l xl Xl + " C l v c L X r XR + v C l v L xr Xl C ipL tpL T 1 " r + h.c. , (40) 



a ip L ip L r r $ + h.c. D -a v L v L X * L xl ®d + av c L v c L X r Xr $d , (41) 
which could come from the Fig. EH and Fig. Eh, together give a neutrino mass matrix: 



where, a = Mp ($£))/m|. This mass matrix is obtained by integrating out the 
heavy modes S^. In the absence of Z)-parity violation, this mass matrix remains 
symmetrical and one of the eigenvalues vanishes, leading to a massless left-handed 
neutrino. When D-parity violating effect is included, the symmetry between the left 
and the right handed neutrinos is lost and the left-handed neutrinos become light 
and massless. In the limit a > — 1 and < irr, diagonalization of this matrix 
gives a light neutrino with mass 



and a similar D-parity violating effective operator 




(42) 




(43) 



m u = 



(1 + a) M s ' 
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This gives the correct order of magnitude for neutrino mass for m L ~ 100 GeV 
and Ms ~ 10 13 GeV. This tiny neutrino mass is of the see-saw type and in fact all 
fermion masses are of the see-saw type in this model. 



r 

® 



R 



r 



(a) (b) 
Figure 3: Tree level diagrams contributing to the neutrino masses. 



5 Gauge coupling unification 

In this section we shall study the renormalization group equations for the evolution 
of the coupling constants in our model. We start with the one-loop renormalization 
group equation for the gauge coupling constants 

where t = ln(/x) where \i stands for the energy-scale of our theory, gi is the gauge 
coupling constant of the group which is a subgroup of the semi-simple gauge 
group G\ x G 2 x • • x G{ x • • • and the beta functions $ contain contributions from 
gauge bosons, fermions and scalars as: 

g 3 

Pi = 1 [{gauge bosons} + {fermions} + {scalars}] . (45) 
To two-loop the (3 functions of any semi-simple gauge group is given as [11]: 

^ a S + T^% b " 9hh (46) 

where the a^s and bijS are the one-loop and two-loop (5 function coefficients respec- 
tively. Here n is the number of groups whose direct product is the semi-simple gauge 
group of the theory, i,j takes on values from 1, 2, • • -n. 
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First we concentrate on the one-loop effect and later we will see the effects of 
the two-loop coefficients on the gauge coupling evolutions. The a^s calculated for 
the various phases are supplied below [11]. 




In the above table Gsm — SU(3) C x SU(2) L x U(l) and the superscripts (SM), 
(LR), (X) indicates the phase in which the numbers are calculated. The renor- 
malization group (RG) equations can now be used to write down the the standard 
model gauge couplings in terms of the SO(10) coupling. Writing ai = gf/4ii, the 
gauge coupling constant matching conditions at the M R scale are: 



1 



a 1Y (M R ) 
1 



- Gsm 



1 



2 

+ 7T- 



Ol2L(M R ) 
1 



a 2L (M R ) 
1 



5a 2R (M R ) 3 (%b-l(M r ) 
1 



G3221 



Gsa 



a 2L (M R ) 
1 

Ol2R{M R ) 

1 

a 3C (M R ) 



G3221 



G3221 



_a 3C (M R ) 

The matching conditions at the M x scale are: 

1 



(47) 
(48) 
(49) 
(50) 



- G3221 



(*b-l{M r ) 
1 

a 2L (M R ) 
1 

a 2 R(M R ) 
1 

a 3C (M R ) 



G3221 



G3221 



G3221 



G3221 



(*4C(M R ) 

1 

a 2L (M R ) 
1 

oc 2 r{M r ) 
1 



G422 



G422 



G422 



G422 



(51) 
(52) 
(53) 
(54) 
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Finally at the Mjj scale, 



1 






1 


_ocw{M v )_ 


G422 




a 10 (Mu)_ 


1 






1 


_a 2L (Mu)_ 


G422 




_a 10 (Mu)_ 


1 






1 


_oi 2 r{M u )_ 


G422 




aio(M [/ )_ 



J SO10 



JSO10 



SO10 



(55) 
(56) 
(57) 



With the help of the above matching conditions and the RG equation we can write 
to one-loop, 



1 



1 



a lY (M z 



a l0 (Mu) 



+ 8tt 



a^y In 



M. 



,(X) 
l 2R 



a 2L (M z ) 



a 10 (M v 

1 



5 

+ 8tt 



a {X) 

a B-L 



hi 



a 3C (M z 



\M X 

+ 87T 



aio(Af, 



!7; 



4f f) ln 



« (5M) In 
a 3C 111 



M z 

Mu 
M x 
Mr 



3 ,(^) , 2 (^)\ 



+ af/)ln 



In 



M 



x 



M, 



M 



x 



Mr 



M 



R 



M, 



+ a 



(LR) 
3C 



hi 



M 



x 



M. 



R 



+ a\ c ' In 



x 



(5* 



(59) 



(60) 



The linear combinations of the gauge couplings that yields sin 2 9w and a s are the 
following: 

1 1 \ 



sin 2 6 W (M Z ) 



1 - 



a(M z 



(61) 

3 a s (M z ) a{Mz) \a 2L (M z ) + 3« 1 F(M Z ) " 3a 3 C(M z )) ' (62) 
where a and a s are related to the electromagnetic and strong interaction coupling 
constants in the present symmetry broken phase. Using the experimental numbers 
[12,13], 

sin 2 9 W (M Z ) = 0.2312, a~ 1 (M z ) = 128.91 and a s (M z ) = 0.1187 . (63) 

Eq. (JoTj) and Eq. (|B12|) reduces to the following: 

1 1 



1 



+ 



a 1Y {M z ) 
5 



«2l(M 2 



29.66, 
106.444 . 



(64) 
(65) 



a 2L (M z ) 3a 1Y (M z ) 3a 3C (M z ) 
The above equations can be utilized for calculating the intermediate scales like Mr 
and Mx in our theory. Here we discuss two cases. 
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When M x = M r 

In this case from Eq. (fHHjl . Eq. (f5T)|) and Eq. and using the (3 function coefficients 
given in the last table we get, 



a 1Y (M z ) a 2L (M z ) 



1 

2^ 



23 

5 \M Z 



5 \M R 



3 In 



29.66. 



(66) 



Similarly from Eq. 
cients we get, 

1 5 



Eq. (O, Eq. and Eq. © and the (3 function coeffi- 



+ 



a 2L (M z ) 3a 1Y (M z ) 3a 3C (M z ) 



1 

2^ 



15 In 

(M v 



Mi 



u 



23 /Mfl 

M z y 3 n VM z 



+ - In . 

3 VM K 



106.444 . 



Eliminating Mu from the above two equations we get 

'Mr 



hi 



M? 



16.301, 



(67) 



(6J 



and if we take M z = 91.1876 GeV then Mr = 1.1 X 10 9 GeV. The above value of M R 
can be taken as the lowest possible value of it in our model and all the predictions 
in our model will be made assuming Mr ^> 10 9 GeV. 



When M x ^ M v 

In this case the two equations corresponding to Eq. 



1 

2^ 



5 \M Z 



1 . 
+ - In 

5 



M 



x 



M, 



R 



5 



Mi 
Mv 



3 In 



and Eq. (jBTj) are: 

M x - 



M? 



29.66 



(69) 



and 



1 

2^ 



3 \M Z 



ilnf^ 

3 \M R 



12 In 



Mr 



M 



x 



15 In 



M 



x 



Mo 



106.444 . (70) 



Eliminating M v from the above two equations gives us a relation between M x and 
Mr as, 



In 



M 



x 



M. 



R 



21 



In 



M, 



R 



M, 



19.214 



(71) 



From the above equation it can be verified that if we take M z = 91.1876 GeV and 
impose Mr = Mx then Mr = 2.01 x 10 10 GeV. In the next subsection we include 
the two loop results and the above results are re-derived computationally. From the 
computational results we see that the above value of Mr ~ 10 11 GeV is two orders 
of magnitude smaller than the actual one. 
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Two-loop result 



After the discussion on gauge coupling unification to one-loop we discuss about the 
two-loop effects of the RG equations. To two-loop the (3 functions are given in 
Eq. (jlBJ). The a^s for the various phases has been supplied in the table appearing in 
the beginning of this section and the b^s for the various phases are as follows: 



ASM) 
ij 



/ 1Z6 
' 25 



11 
10 



18 
5 



Here i,j = 1Y, 2L, 3C, 



u(LR) 
ij 



( 

3 
3 
VI 



8 9 





9 
2 



9 




and here i,j 



B — L } 2L } 2R } 3C, 



/ 199 
' 3 



V ¥ 



199 

3 

21 
2 



44 \ 
5 » 



12 



-26 J 



2 

12 
12 
-26 / 



105 
2 

105 
2 



1307 , 
6 / 



(72) 



(73) 



(74) 



Where i,j = 2L, 2R, AC. 

If we start from jj = Mj< = 173 GeV (where n f = 6) and fix M R = 10 11 GeV 
then the evolution of the gauge coupling constants are as given in Fig. 0] In the next 
phases the coupling constant evolution shows that at Mx both a^c and oib-l unite 
to produce a^c- Fig. |U shows that from Mx on wards the development of a.2L and 
CK2R are identical. At around M\j ~ 10 15 GeV the gauge coupling constants unite. 
Our computational results show that the highest value of Mr is just slightly above 
10 11 GeV. If Mr is much above the above mentioned value then Mx comes down 
and M R >M X . 

In models with triplet Higgs scalars, it is possible to consider Mr > Mx- 
At the scale Mr, the group SU(2) H is broken into 11(1)^ and stays orthogonal to 
SU(4) C . Later when SXJ(4:) C is broken to 811(3)^ x U(1) B _ L . Subsequently at a 
much lower scale the symmetry breaking, 11(1)^ x U(1) B _ L — ► U(l) y , takes place. 
However, this is not possible in scenarios with only doublet Higgs scalars, since the 
right-handed Higgs scalar doublet does not have any component with B — L quantum 
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Figure 4: Plot showing the evolution of a^ 1 s in the SM phase, left-right phase upto M\j. 
The abscissa is ln(/z), where initial is M t = 173 GeV and Mr = 10 11 GeV in our case. The 
figure in the right shows explicitly the gauge coupling unification at Mjj ~ 10 15 GeV. 



number to be zero and hence it breaks SU(2) H and \]{1)b-l simultaneously. As a 
result, the highest value of M R could be 10 11 GeV. As we discussed earlier [14], the 
gauge coupling unification also requires a lower bound on Mr > 10 9 GeV. From this 
lower (higher) bound for Mr, we estimated computationally higher (lower) bound 
for Mx and Mjj as shown in Figs. El 13 The lower and higher bounds for Mx 
are ~ 10 11 GeV and ~ 5 x 10 11 GeV respectively. Also lower and higher bounds 
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for Mu are estimated as 5 x 10 14 GeV and ~ 10 15 GeV respectively. In this regard 
it is important to note that the stability of the proton offers further constraints on 
the parameter space of their model, since the GUT scale is lower than most of the 
conventional models and is smaller than 10*- 15) GeV. 

The difference between Mj and Mr can go up from (i.e, Mr = Mx\ no 
intermediate G3221 symmetry) to ~ 550 GeV as shown in Fig. |S| This range of 
Mr plays an important role for Yukawa coupling unification and leptogenesis in this 
model. 



6 Yukawa coupling evolution 

In sec. IH Eq. ffTTTTf) gives the only Yukawa couplings of our theory. Eq. (J33|) is valid 
in the SO(10) level. In the G422 level the Yukawa couplings will become: 

4 4 * 2 ' 2) = g'S L FP^ L + g"S L Fi r) ^ R + M S S L S L + h.c. . (75) 
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Figure 6: M R vs. My shows lowest and highest bound. 



The gauge transformation properties of the various fields, except Sl which is a gauge 
singlet, are shown below: 



F® = (4,2,1) C 16 = T, 

FP EE (4,1,2) C 16 EE T, 

$2 EE (4,1,2) C 16 EE Tt, 

$ fl EE (4,1,2) C 16 EE rt. 



(76) 



In the above equations (x, y, z) designates the transformation properties under G422. 
The Yukawa couplings in the G3221 phase is as: 



4 3 * 2 ' 2,1) = k[ q) S L q L Hl + k' {l) S L £ L xl + k'l q) S L {q R ) c H R + k^S L (i R ) c x R 
+ M s S L S L + h.c. 



(77) 



The gauge transformation properties of qi, £l, (qr) c , (^r) c , Xl an d Xr are specified in 
Eq. (jlU|) , Eq. (|T2*|) and Eq. Q . Here we specify the gauge transformation properties 
of the other two fields present in the above equation. 



HI 

Hr 



(3,2,1,-i) C (4,2,1) C 16 ee Tt, 
(3,1,2,|)C(4,1,2)C16 ee T. 



(78) 



In the above equations (x, y,z, w) designates the G3221 transformation properties 
and (x, y,z) designates the G422 transformation properties. 
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In the next stage, that is in the G321 phase, the Yukawa couplings look like: 

£y 2,1) = V^SlQlTI + y{iL)S L L L h* L + V{u)S l {Qr) c u T r ~> 
+ y {d) S L {Q R ) d T R +) + y {e) S L (erf <j><$ 
+ y{u)S L {v R ) c ^ ) R + M S S L S L + h.c. . (79) 

Here the various standard model fermion fields transform under G321 as: 

Ql = (3,2,|)C(3,2,1,|)C(4,2,1)C16 = T, 

L L = (1 2,-1) C^l,2, 1,-1) C(_4, 2,1) C 16 = T, 

(Qr)u = (3,1,|)C(3,1,2,-|)C(4,1,2)C16 = I\ 

(Q R y d = (3,1,-|)C(3,1,2,-|)C(4,1,2)C16 = T, { ™> 

(e R ) c = (1,1,-2) C (1,1,2,1) C_(4, 2,1) C 16 = T, 

(is R ) c = (1,1,0)C(1,1,2,1)C(4,2,1)C16 = r. 

Similarly the various Higgs fields transform as: 



Ti = 


(3,2,-|)c(3,2,l,-|)c(4,2,l)cl6 


= rt, 


h* L = 


(1,2,1)C(1,2,1,1)C(4,2,1)C16 


= rt, 




(3,1,-|)C(3,1,2,|)C(4,1,2)C16 


= rt, 




(3,1,|)C(3,1,2,|)C(4,1,2)C16 


= rt, 




(1,1,2)C(1,1,2,-1)C(4,2,1)C16 


= rt, 




(1,1,0)C(1,1,2,-1)C(4,2,1)C16 


= rt. 



(81) 
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Figure 8: Mx — Mr vs. Mx shows how Mx — Mjj approaching zero 



In the above expressions the first triplet (x, y,z) designates the transformation prop- 
erties under G 32 i, the next four numbers (x, y, z, w) designates the G3221 transforma- 
tion properties and the triplet (x, y, z) stands for the G422 transformation properties. 

Now we give an order of magnitude estimation about the running of the effec- 
tive top Yukawa coupling in our theory. From Eq. (j3*U|) it is seen that the effective 
top quark Yukawa coupling looks like (X' F v rP) / '(87r 2 ] M X ) ■ In the G 32 i phase it looks 
like {\' T VRy^ qL ^) / {9>ti 2 Mx) in the convention adopted to name the Yukawa couplings 
in Eq. ()79|) . Calling this effective coupling as y t it will evolve simply like: 



dyt 
dt 



Vt 



(16vr 2 



(82) 



as in the standard-model up to the Mr scale. The gauge and quartic coupling con- 
tributions will be negligible compared to y t . Starting from the top-quark mass at 
the electroweak scale, the evolution equation gives the effective top-quark Yukawa 
coupling at the left-right symmetry breaking scale to be of the order of 0.76 for 
Mr ~ 10 11 GeV. Since the effective coupling constant is a product of three couplings 
(^r v RV(qL)) I '(87r 2 Mx) and if vr ~ M x then each of these couplings can individu- 
ally take values large enough as 3.9. As a result the Yukawa sector becomes non- 
perturbative in the G3221 phase. But on the other hand if we have (Mx/vr) ~ 10~ 2 
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then the situation changes. In this case the individual coupling becomes of the order 
of 0.84 which implies the theory is still perturbative. Mx can be the heavy gauge 
boson masses or the singlet fermion Sl mass. The previous condition means that 
the heavy gauge bosons or the singlet fermion cannot be heavier than 10 9 GeV in 
our theory if we take Mr ~ 10 11 GeV for a perturbative scenario in the Yukawa 
sector. 

Above the left-right symmetry breaking scale Mr up to the unification scale 
Mu, the coupling constants A' r and / will evolve separately. The separate Yukawa 
couplings remains finite up to the unification scale. 

7 Leptogenesis 

Since the neutrino masses now depend on the couplings with the singlets, there is 
no stringent restriction coming from the up quark masses. As a result, it may be 
possible to get large neutrino mixing angles. The right-handed neutrinos and the 
new singlet fermions can now decay into light leptons. The Majorana masses of 
the left-handed and right-handed singlets violate lepton numbers, which in turn can 
generate enough lepton asymmetry. Before the electroweak phase transition this 
asymmetry can then generate a baryon asymmetry of the universe [15]. Since there 
is no supersymmetry, the gravitino bounds are not present. The out-of-equilibrium 
condition can be satisfied near the GUT scale since the couplings are large to get 
the required neutrino mass with large see-saw scale. In this model there is another 
interesting feature that the singlets combine with the right-handed neutrinos to form 
pseudo-Dirac particles and hence resonant leptogenesis may also be possible [16,17]. 

For leptogenesis, consider the interactions of Eq. (j3HJ). Unlike usual see-saw 
models with triplet Higgs scalars [18, 19], in this model the right-handed neutrinos 
cannot decay into left-handed neutrinos and Higgs bi-doublets dominantly. The 
simplest lepton number violating interactions come from the decays of Sf 

Si -> £jL + X*L> 

-> ijL+XL- (83) 

The Majorana masses of Si allow the singlet to decay into both leptons and antilep- 
tons violating lepton number by two units. In the present model both £l and xl 
are very light and hence these decays are allowed. 

For CP violation there are two types of one-loop diagrams which interferes with 
the tree-level diagrams for the decays of S^ These are the vertex type diagrams of 
Fig. El and Fig. EH 
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The right-handed neutrinos do not take part in leptogenesis directly, but due 
to mixing of the right-handed neutrinos with the heavy singlets Si, the right-handed 
neutrino decays also enter in the picture of leptogenesis. In the limit of M s ^> m R , 
the right-handed neutrinos and the Si are both heavy and distinct. In this case the 
amount lepton asymmetry due to CP violation is given by, 

1 MiM£ Q (fc(o al^(0a2) 2/3(^(0 /31^(Z) /32)] 



8ttM 2 EaWn 



12 



(84) 



->a I ^(7) all 

where we assumed that the mass matrix of 5, are diagonal and the eigenvalues are 
hierarchical M 1 < M 2 < M 3 where Mj are masses of We can write the effective 
lepton asymmetry as: 

e = 5(, (85) 
where ( is a suppression factor which depends on the amount of departure from 




Figure 9: Vertex type diagrams interfering with tree level diagram. This is similar 
to the CP violation coming from the penguin diagram in K-decays. 




Figure 10: Self energy diagram interfering with tree level diagram. This is similar 
to CP violation in K — K oscillation, entering in mass matrix of S^. 
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equilibrium. The exact value of ( can be obtained by solving the Boltzmann equation 
taking all the interactions into consideration. However, it is also possible to make an 
order of magnitude estimate for the amount of asymmetry, which will be very close 
to the actual value. Since Si is the lightest of the singlets, the decay of this singlet 
will be able to generate the lepton asymmetry. The asymmetry generated or washed 
out by the heavier ones S 2 or S3 will be smeared out by the interactions of Si after 
S 2 and S 3 had decayed away. So, for an estimate we shall only consider the decays 
of Si. This assumption is well justified when the singlets Sj have a hierarchical mass 
structure. 

The out-of-equilibrium condition is parametrized by: 



When 9 > 1, there is no Boltzmann suppression of the generated asymmetry and 
the out-of-equilibrium condition is satisfied. In this case the generated asymmetry 
is given by 5 and it is not washed out after it is created and one gets ( = 1. 
However, if 9 3> 1, then the interaction strength is so slow that the generated 
asymmetry can never reach the value 5. Although the interactions cannot wash out 
the asymmetry after it is generated, the amount of asymmetry is less than S. In the 
case of 9 < 1, the generated asymmetry is same as the CP asymmetry of 5, but even 
after the asymmetry is created, the interaction remains strong enough to deplete the 
asymmetry. Although the depletion is exponentially fast, it cannot compete with 
the expansion of the universe for long and the final amount of asymmetry is not 
exponentially depleted. It was shown that [20] the suppression factor ( is almost 
linearly proportional to 9. In the present model we come across this last scenario. 

In the present model 10 9 GeV < M R < 10 11 GeV. While a lower value of M R 
is preferable for out-of-equilibrium condition, since the Yukawa couplings grow very 
fast above the scale M R we have to consider the highest value of M R . Taking the 
hierarchical structure of Sj, we consider the mass of the lightest singlet to be around 
10 10 GeV. Taking M 1 ~ 10 10 GeV, g* ~ 10 2 and k[ l)al ~ 1 we find that 9 is much 
lower than 1, which gives a strong suppression factor of ( ~ 1CT 7 . On the other 
hand the Yukawa couplings in this model comes out to be of the order of 1 and 
hence we get a large enhancement in the CP asymmetry and 5 in our case can be 
as large as 10~ 2 and so the lepton asymmetry parameter e ~ 1CT 9 . 

At this stage, AB = and AL is given by 5. Thus 




(86) 



A(B-L) = e. 



(87) 
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The final baryon asymmetry after the electroweak phase transition is thus given by 

1^ = I A (5 - L) = - ~ 1(T 10 . (88) 
s 3 3 

In the case of Ms <C m,R, the right-handed neutrinos and the Si singlets of 

every generation are almost degenerate. The mass splitting between the states vm 

and Si with mass m R is of the order of M s . Although Si decays will now generate a 

lepton asymmetry, both the heavy mass eigenstates contain the states Si. As a result, 

when these two almost degenerate states decay, there may be resonant leptogenesis 

(which will however require new interactions and fine-tuning) and hence the scale 

of leptogenesis could be very low. For the present scenario since Mr and hence Mi 

cannot be much lower, this is not important and hence we shall not discuss it in any 

further detail. 

8 Summary 

In conclusion, we constructed an SO(10) GUT without any Higgs bi-doublets. All 
the symmetry breaking could be achieved by only two Higgs scalars, a 210 and a 
16. By including a massive singlet fermion per generation we break chiral symmetry 
which can then give masses to all the fermions radiatively without introducing any 
new scalar fields. All fermion masses have the same see-saw form. The spontaneous 
parity breaking plays a crucial role in breaking the left-handed and right-handed 
SU(2) groups at two widely different scales and also giving masses to the left-handed 
neutrinos in this scenario. The spontaneous breaking of an ungauged discrete sym- 
metry, the D-parity, which is a special feature of this model may cause formation 
of very heavy domain walls of GUT-scale mass. Since the GUT scale in the model 
is rather low (~ 10 15 GeV), some of the parameters of the model may have to be 
constrained further to prevent fast proton decay. These features require further 
study and will be taken up in future. The model allows large neutrino mixing and 
required neutrino masses. The baryon asymmetry of the universe can be explained 
through leptogenesis. 
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